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ABSTRACT 

Context. This paper extends the models of Craig & McClymont (1991) and McLaughlin & Hood (2004) to include finite fi and 
nonlinear effects. 

Aims. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. 

Methods. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 
2001) and consider an initial condition in v x B • z (a natural variable of the system). 

Results. We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a 
'cusp-like' point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current 
sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but 
in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of 
the final current density left in the system. 

Conclusions. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime. 
Key words. Magnetohydrodynamics (MHD) - Waves - Shock waves - Sun: corona - Sun: magnetic fields - Sun: oscillations 



1. Introduction 

, It is now known that MHD wave motions (e.g. Roberts [20041 De 

■ Moortel|2003 Nakariakov & Verwichte [2()(BI are omnipresent 

■ throughout the solar corona (Tomczyk et al. 120071 ). Many so- 
lar instruments have observed various MHD wave motions in 
the solar atmosphere: slow magnetoacoustic (MA) waves have 
been seen in SOHO data (Berghmans & Clette |1999. Kliem 

! et al. 120021 Wang et al. |2U02] l and TRACE data (De Moortel 
' et al. I2000I I. Fast MA waves have been seen with TRACE 
; (Aschwand en et a l. [19991 [20021 Nakariakov et al. [1999 1 Wang 
& Solanki [2004l l and Hinode (Ofman & Wang 2008). Non- 
thermal line narrowing / broadening due to Alfven waves has 
been reported by Harrison et al. (I2002I I / Erdelyi et al. ( |19981 l 
and O'Shea et al. (2003) . More recently, Alfven waves have 
been observed in the corona (Okamoto et al. 120071 Tomczyk et 
al. 120071 ) and chromosphere (De Pontieu et al. 120071 ). although 
these claims are currently subject to intense discussion (Erdelyi 
& Fedun 2007; Van Doorsselaere et al. [2008l ). 
It is clear that the coronal magnetic field plays a fundamental 
role in the propagation and properties of MHD waves, and to be- 
gin to understand this inhomogeneous, magnetised environment, 
it is useful to look at the topology (structure) of the magnetic 
field itself. Potential-field extrapolations of the coronal magnetic 
field can be made from photospheric magnetograms, and such 
extrapolations show the existence of important features of the 
topology: null points - locations in the field where the magnetic 
field, and hence the Alfven speed, is zero, and separatrices - 
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topological features that separate regions of different magnetic 
flux connectivity. Detailed investigations of the coronal mag- 
netic field, using such potential field calculations, can be found 
in e.g. Brown & Priest (2001 ), Beveridge et al. (2002), Regnieret 
al. ( 120081 or a more comprehensive review by Longcope ( 120051 ) . 

The propagation of fast magnetoacoustic waves in an inhomo- 
geneous coronal plasma has been investigated by Nakariakov & 
Roberts ( 1995 ), who showed that the waves are refracted into re- 
gions of low Alfven speed. In the case of null points, the Alfven 
speed actually drops to zero. 

McLaughlin & Hood (2004) (hereafter referred to as [Paper Ij ) 
solved the linearised, = MHD equations using a two-step 
Lax-Wendroff numerical scheme. They found that in the neigh- 
bourhood of a single 2D X-point, the fast MA wave refracted 
around and accumulated at the null point. These key results have 
been found to carry over from the simple 2D single magnetic 
null point to two null points (McLaughlin & Hood 2005 ) and to 
a more realistic magnetic configuration of a null point created 
by two dipoles (McLaughlin & Hood 2006al ). It should be noted 
that the behaviour of the fast wave is entirely dominated by the 
Alfven-speed profile, and since the magnetic field drops to zero 
at the null point, the wave will never reach the actual null for a 
/3 - plasma. McLaughlin & Hood ( 2006b ) extended the model 
of [Paper I to include plasma pressure effects. This naturally led 
to the inclusion of the slow MA wave and the introduction of a 
y6 = 1 layer around the null point; representing a high /3 environ- 
ment inside and low outside. Coupling and mode conversion is 
observed at locations where the sound speed and Alfven speed 
become comparable in magnitude (e.g. Zhugzhda & Dzhalilov 
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[19821 Callv 120011 Bloomfield et al. l2006l McDougaU & Hood 
I2007l l. 

Waves in the neighbourhood of a single 2D null point have 
been investigated by various authors. Bulanov & Syrovatskii 
( 1198 0^ provided a detailed discussion of the propagation of har- 
monic fast and Alfven waves using cylindrical symmetry. Craig 
& Watson (1992) mainly considered the radial propagation of 
the m -Q mode (where m is the azimuthal wavenumber) using a 
mixture of analytical and numerical solutions. They showed that 
the propagation of the m - wave towards the null point gen- 
erates an exponentially large increase in the current density and 
that magnetic resistivity dissipates this current in a time related 
to log rj. Craig & M cClymont ( 1 99niT993T l. Hassam (^19921) and 
Ofman et al. ( 119931 ) investigated the normal mode solutions for 
both m - and m + modes with resistivity included. Again, 
they emphasise that the current builds up as the inverse square 
of the radial distance from the null point. All these investigations 
were carried out using cylindrical models in which the generated 
waves encircled the null point. 

Reconnection can occur when strong currents cause the mag- 
netic fieldlines to diffuse through the plasma and change their 
connectivity (Parker [19571 Sweet [19581 Petschek [T964l l. In 2D, 
reconnection can only occur at null points (Priest & Forbes 
12000) . Dungey (119531) reported that a perturbed X-point can col- 
lapse if the footpoints of the field are free to move, Mellor et 
al. ( 12002 ) looked at the linear collapse of a 2D null point, and 
Imshennik & Syrovatsky ( 119671 ) described the collapse with an 
exact, non-linear solution of the ideal MHD equations. However, 
these papers did not include the effect of gas pressure, which 
would act to limit the growth of the current density. In con- 
sidering the relaxation of a 2D X-type neutral point disturbed 
from equilibrium, Craig & McClymont (1991) found that free 
magnetic energy is dissipated by oscillatory reconnection, which 
couples resistive diffusion at the null to global advection of the 
outer field. An example of oscillatory reconnection generated by 
flux emergence within a coronal hole was recently detailed by 
Murray et al. (2008). Finally, Longcope & Priest (2007) inves- 
tigated the diffusion of a 2D current sheet subject to suddenly 
enhanced resistivity. They found that the diffusion couples to a 
fast MA mode which propagates the current away at the local 
Alfven speed. 



The aim of this paper is to extend the model used in Paper I 
and Craig & McClymont ( 119911 ) to include finite (3 and nonUn- 
ear effects. To realise this, we will be solving the compressible 
and resistive MHD equations using a Lagrangian remap, shock 
capturing code: LARE2D (Arber et al. 120011 ). The key results 
from [Paper I| i.e. for the linear fast wave, are that it demon- 
strates refraction, that the wave energy accumulates at the null, 
and that current density builds up exponentially at that point. We 
believe it is important to extend this work to include nonlinear 



effects since Paper I indicated a preferential topological location 



for (ohmic) heating and we need to see if this observational pre- 
diction persists when we consider larger (and hence nonlinear) 
wave amplitudes. This paper will address three main questions 
that naturally arise when extending faper I| into the nonlinear 
regime: 

(1) Does the fast wave now steepen to form shocks, and can 
these propagate across or escape the null? 

(2) Can the refraction effect drag enough magnetic field into the 
null to initiate X-point collapse or reconnection? 

(3) Has the rate of current density accumulation changed, and is 
the null still the preferential location of wave heating? 




Fig. 1. The equilibrium magnetic field. The red lines denote the 
separatrices and arrows indicate the direction of the magnetic 
field. The null point is located at the origin, where the separatri- 
ces intersect. 



The paper has the following outline: the basic setup, equations 
and assumptions are described in ^ ^details the nonlinear fast 
MA wave behaviour and the resultant oscillatory reconnection is 
discussed in ^ We consider different amplitudes for our initial 
condition in ^and the conclusions are given in ^ 



2. Basic Equations 

We consider the 2D compressible and resistive MHD equations 
appropriate to the solar corona: 



P 



— -^(vV)v 
of 

5B 



= -VpH-l-VxBlxB , 



— = V X (v X B) -H T/V-'B 
ot 



dp 

^+V-(pv) = 0, 
ot 



dt 



+ {y-V)e 



(T 



(1) 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field), p is the 
plasma pressure, ^ - Anx lO^^Hm"' is the magnetic permeabil- 
ity, cr is the electrical conductivity, rj - \ /ficr is the magnetic 
diffusivity, e = p/piy- I) is the specific internal energy den- 
sity, where y = 5/3 is the ratio of specific heats and j = V x B/yU 
is the electric current density. 

The LARE2D numerical code utilises artificial shock viscosity to 
introduce dissipation at steep gradients. The details of this tech- 
nique, often called Wilkins viscosity, can be found in WiUcins 
( [T9801 ) and Arber et al. ( i200Tl ). 
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2.1. Basic equilibrium and non-dimensionalisation 

The equilibrium magnetic field structure is taken as a simple 2D 
X-type neutral point. The aim of studying waves in a 2D config- 
uration is one of simplicity. There are many complicated effects 
including mode conversion and coupling, and a 2D geometry al- 
lows us to understand and explain these behaviours better, before 
the extension to 3D. The initial magnetic field is taken as 



Bo = ^(3',^,0) , 



(2) 



where B is a characteristic field strength and L is the length scale 
for magnetic field variations. This magnetic field can be seen in 
Figure[T] Equation (|2]i is slightly different to that used in |Paper l"} 
where the authors considered Bq - (x, 0, -z). This simply rep- 
resents a ;r/4 rotation of our magnetic field and the key results 
are still valid. Note that this magnetic configuration is no longer 
valid far from the null point, as the field strength tends to infin- 
ity. However, McLaughlin & Hood (I2006al l looked at a yS = 
magnetic field that decays far from the null and they found that 
the key results from |Paper I| remain true close to the null. 
We will also find it useful to consider A, the vector potential, 
where B = V x A. In 2D and for the coordinate system used in 
this paper, A = A^z and our equilibrium vector potential is given 
by: 



Ao = A()Z = -x^)z 



(3) 



We now consider a change of scale to non-dimensionalise all 
variables. Let v = vqv*, B = BB*, x - Lx*, y = Ly*, p = pop*, 
P - PoP*7 V = 7^*^ t - tot*, A = BLA* and rj - jjo, where 
we let * denote a dimensionless quantity and vo, B, L, po, po, to 
and 7/0 are constants with the dimensions of the variable they are 
scahng. We then set B/ -sjppo - vq and vq - L/tg (this sets vq 
as a constant background Alfven speed). We also set rjotolL^ - 
R^^, where R„, is the magnetic Reynolds number, and set Po - 
2fipo/B^, where /3o is the plasma-j6 at a radius L from the origin. 
This process non-dimensionalises equations ([U and under these 
scalings, f* = 1 (for example) refers to f = fo = L/vq; i.e. the 
time taken to travel a distance L at the background Alfven speed. 
For the rest of this paper, we drop the star indices; the fact that 
all variables are now non-dimensionalised is understood. 
We take the equilibrium density to be uniform, i.e. p - po; 
a spatial variation in po can cause phase mixing (Heyvaerts & 
Priest , 1983. De Moortel et al. .l999i Hood et aL l2002] l. We set 
Rm = 10^. 

Finally, we consider the equilibrium plasma to be cold: T - OK 
(i.e. jSo = 0) and, hence, ignore plasma pressure effects (as in 
[Paper However, as we will see below, magnetic shocks heat 
the plasma and so the plasma will not remain cold (see e.g. §1.5 
in Priest & Forbes 120001 1 . 



from the null point and, hence, waves accelerate as they propa- 
gate outwards. Since we do not want reflected waves to influence 
our null point, implementation of such a damping region is es- 
sential. 

As seen in [Paper H there are two natural variables to consider in 
our system: Vj_ = (v x B) ■ z = Vi-Z?,. - VyB.t and vy = v ■ B = 
VjB.v + VyBy. Here, v^ and vy are related to the perpendicular and 
parallel velocity, respectively and, as seen in Paper I , their im- 
plementation naturally simplifies the governing equations, aids 
in MHD mode interpretation and (for Vj^) led to an analytical 
solution to the linear, cold plasma equations. 
In cartesian coordinates: v y = (vyBy + v^By) / |B|^ and Vy = 
(v||By - Vj^By) / |B|^. In polar coordinates: Vy = v^ cos 6/r, Vy = 

-Vj_ sin 6/ r, where r - ^Jx^ + y^ and we take vy to be initially 
zero. We note that with our choice of magnetic null point (equa- 
tion|2]), if we drive any of the velocity variables, the system will 
naturally develop a dependence. 

Paper I clearly demonstrates that the Alfven speed (v^ = B^ -n 
By = -H y^) plays a vital role. Hence, it is natural to consider 
either a polar coordinate system or a circular pulse. In addition, 
as commented by McClements et al. ( I2004I I. a disturbance ini- 
tially consisting of a plane wave is refracted as it approaches the 
null in such a way that it becomes more azimuthally symmetric. 
Thus, it is appropriate to consider the evolution of azimuthally 
symmetric perturbations. We expect the nonlinear behaviour to 
be more complicated than the linear equivalent and so, in order 
to clearly demonstrate the differences between the two systems, 
in this paper we consider an initial condition in velocity, such 
that: 

v^(x,y,t = 0) ^ 2Csm[7T(r- 4.5)] for 4.5 < r < 5.5 , (4) 
vy(jc,y,f = 0) = 

where 2C is our initial amplitude. Initial condition (|4| describes 
a circular, sinusoidal pulse in v^. When the simulation begins, 
this initial pulse will naturally split into two waves, each of am- 
plitude C, travelling in different directions: an outgoing wave 
and an incoming wave. In this paper we will focus on the incom- 
ing wave, i.e. the wave travelling towards the null point. The 
damping regions will remove kinetic energy from the outgoing 
waves, and so they do not influence the null. The initial condi- 
tion produces a propagating disturbance that crosses magnetic 
fieldlines. Hence, we identify these waves as (initially) fast MA 
waves. 

By choosing a small value for C in equation (|4|, we can recover 
the linear results from [Paper I This can be seen in Appendix [A] 
where we set C = 0.001. For the nonlinear work described in ^ 
and ^ we set C = 1 . Different values of C will be considered 



2.2. Initiai and boundary conditions 

Equations ^ are solved numerically using a Lagrangian remap, 
shock-capturing code called LARE2D (Arber et al. 2001). The 
equations are solved computationally in a square domain x,y e 
[-20, 20] with a numerical resolution of 5120 x 5120. Zero gra- 
dient boundary conditions are applied to the variables B, p, e at 
the four boundaries, and v is set to zero on all boundaries, i.e. 
reflective boundaries. A damping region exists for x^ + y^ > 6 
and so all oscillations that enter this region are slowly damped 
away. The (equilibrium) Alfven speed increases with distance 
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Fig. 2. Contours of Vj. for a fast wave pulse initially located at a radius r - 5, and its resultant propagation at (Alfven) times 
f = 0,0.1,0.5 & 1. The black lines denote the separatrices and the null point is located at their intersection (the origin). The full 
evolution, < f < 60, is available as an mpeg animation in the online edition of the Astropnomy & Astrophysics Journal. 
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Fig. 3. (a) Plot of Vi.(x,0) for -5 < x < 5 at f = 1. (b) Plot 
of Vv(0,3') for -5<3'<5atf=l. The red line indicates the 
location of the null point. For x,y > 0, we have plotted stars to 
indicate the grid resolution. 



3. Nonlinear fast MA behaviour 

3.1. Development of shocks (0 < f < 1) 

The evolution of Vj^ can be seen in Figure |2](0 < t < 1), Figure 
IS (1.4 < / < 2.8) and Figure |6] (2.9 < f < 60), and the fuU 
evolution, < f < 60, is available as an mpeg animation in the 
online edition of the Journal. 

From Figure |2] we see that the initial pulse has split into two 
oppositely travelling wave pulses, where both waves can be seen 
at f = 0.1. In order to clearly show the wave behaviour, we have 
plotted x,y e [-5,5], i.e. a square subsection of our total com- 
putational box. Hence, only the incoming wave can be seen af- 
ter time t = 0.1. This figure can be compared with Figure [ATI 
which corresponds to the linear regime. 

There are two key features to note from Figure|2] Firstly, we see 
that the incoming wave propagates across the magnetic fieldlines 
and keeps its initial pulse profile, i.e. an annulus, and the maxi- 
mum amplitude remains at C = 1 . The annulus contracts as the 
wave approaches the null point and this is the same refraction ef- 
fect described in 



Paper II This refraction effect occurs since the 



(equilibrium) Alfven speed is spatially varying. 
Secondly, we note that the incoming wave pulse is develop- 
ing an asymmetry: the wave peaks are propagating faster (rel- 
ative to the footpoints) in the y-direction than the x-direction. 
Thus, the wave pulse is developing discontinuities, where in the 
y-direction the wave peak is catching up with the leading foot- 
point, and in the x-direction the trailing footpoint is catching 
up with the wave peak. These discontinuities can be clearly seen 
in Figure (trailing edge in Vj^[x, 0]) and Figure |3j) (leading 



edge in Vv[0,3']). We have also indicated the grid resolution (us- 
ing stars) on the right-hand side of each subfigure which shows 
that the developing discontinuity is well resolved. 
This development of discontinuities was not reported in ' Paper I| 
as it is an entirely nonlinear effect, which arises because of our 
choice of a velocity initial condition (equation HJ. In the non- 
linear regime, specifying an initial velocity condition also pre- 
scribes a background velocity profile, and this profile can be 
seen in Figure |4] It is this background velocity profile that leads 
to the development of discontinuities on the leading edges in 
the 37-direction and on the trailing edges in the x-direction. 
The Vf (red arrows) and v,, (blue arrows) contributions to this 
background velocity profile can be seen in Figure |4j;. This phe- 
nomenon is most easily understood by means of a simple ID 
example, details of which can be found in Appendix iB] Note 
that the background velocity profile prescribed by equation (HJi 
appears as the m - mode in v^ but corresponds to the m - 2 
mode in cartesian components. 
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Fig. 4. Our choice of initial condition for (a) v v and (b) Vy at f = 0. This choice of initial condition prescribes a background velocity 
profile, (c) Cartoon representation of the Vi (red arrows) and v,, (blue) contributions to this background velocity profile, where the 
black lines denote the separatrices and the null point is located at their intersection. 



3.2. X-point collapse (1.4 < f < 4) 

The wave evolution between 1 .4 < f < 2.8 can be seen in Figure 
|5] From the first row, 1.4 < f < 1.8, we see that the asymmetry 
seen in Figure|2]has now lead to the formation of shock waves. 
From the second and third rows (2 < f < 2.5) of Figure|5] we see 
that the shocks above and below y - have started to overlap 
(starting at x « ±1). This overlap leads to the development of 
hot jets. 

The white box in Figure |5] at f - 2 is analysed in Figure [T^. 
By considering the physical quantities along a line perpendicu- 
lar to the shock front, we see that there is an abrupt increase in 
density, temperature and, consequently, pressure (FigurelTJ?). 
increases in magnitude, whereas By decreases, and both preserve 
their original directions (Figure |7};). Hence, the shock makes 
B refract away from the normal and so we identify it as a fast 
oblique magnetic shock. In the idealised limit v || B, this would 
be called a switch-on shock (see e.g. §1.5.2 in Priest & Forbes 
I2000I I. It is interesting to note that the shock has heated the 
plasma and so y6 in these locations. 

The white box in Figure |5] at f = 2.4 is analysed in Figure [S] 
Figure[8^ shows a contour of v^ for 0<;c<2, -1 < y < 1, and 
we see that the overlap of the shock waves forms a triangular 
'cusp', which we call the shock-cusp. Figure[8]3 shows that these 
hot jets heat the plasma, substantially more than the previous 
shock heating (compare magnitudes of T from Figure |7}3 and 
Figure [S}?). These hot jets also significantly bend the magnetic 
fieldlines. 

The hot jets set up new shock waves emanating from the shock- 
cusp, and by considering the changes of physical quantities per- 
pendicular to the shock front (similar to before), we see that there 
is an abrupt increase in B^ and decrease in By, although both 
preserve their original directions (Figure [8]:). Thus, the shock 
makes B refract towards the normal, and we identify this as a 
slow oblique magnetic shock. In the idealised limit v || B, this 
would be called a switch-off shock. 

The structure of the hot jet is in good agreement with that de- 
scribed by Forbes ( 1988) (see his Figure 15). In addition to the 
slow shocks downstream of the tip of the jet, there is evidence 
of slow shocks along the sides of the jet upstream of the tip (see 
Figure [8J)). There are also kinks in the fieldlines at the tip of 
the jet, indicative of a fast oblique magnetic shock. Thus, the jet 
heating itself is accomplished by a combination of slow and fast 
shocks. The abrupt jump in temperature and magnetic field at 



the tip of the jet is also consistent with the jump predicted by 
Soward and Priest (fT982] l. 

The fourth row of Figure |5] shows that the shocks, both the fast 
shocks and their overlap (the tails of the jets) have reached the 
null point (at t = 2.6). The shocks have deformed the magnetic 
field such that the separatrices now touch one another rather than 
intersecting at a non-zero angle (called 'cusp-like' by Priest & 
Cowlev lI975l l. That the perturbation can reach the neutral point 
is a phenomenon not seen in Paper I (where the perturbation 
reached the null after an infinite amount of time). 
After time t - 2.6, we see that the shock wave can now pass 
through the null point: entirely different behaviour to that seen 
in |Paper I ' This is more clearly seen in Figure |9^, which shows 
a plot of v^,(0,3') at f = 2.6 (just before null point is crossed) and 
t - 3.0 (after null has been crossed). 



6 



McLaughlin et al.: Nonlinear fast wave propagation around an X-point 




Fig. 5. Contours of at various (Alfven) times between t - 1.4 and t - 2.8. Note that the amplitude varies substantially throughout 
the evolution, and hence each row is assigned its own colour bar The black lines denote the separatrices. The white boxes at f = 2 
& 2.4 are analysed in Figures |2]& [8] respectively. 




Fig. 6. Contours of at various (Alfven) times between t - 2.9 and t - 60. Note that the amplitude varies substantially throughout 
the evolution, and hence each row is assigned its own colour bar. The black lines denote the separatrices. The red lines at f = 60 
denote the potential separatrices. 
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Fig. 7. (a) Contour of at f = 2 for 0.2 < x,^ < 0.5. The black 
line denotes the separatrix, the dotted, y - 4.382x - 0.938, and 
dashed, y = -0.228x + 0.341, white lines denote the lines per- 
pendicular and parallel to the shock front, respectively, (b) Plots 
of Temperature, T (red). Pressure, P (blue), and density, p/100 
(green), and (c) plots of B-c (red) and By (blue) perpendicular to 
the shock front. 



Fig. 8. (a) Contour of v_l at f = 2.4 for < x < 2, -1 < 
y < 1. The black line denotes the separatrices, the dotted, 
y = -l.lllx + 1.427, and dashed, y = 0.9x - 0.685, white lines 
denote the lines perpendicular and parallel to the shock front, re- 
spectively, (b) Contour of temperature at f = 2.4 for < x < 2, 
-1 < y < 1. The white lines denote magnetic fieldlines (thick 
white line denotes the separatrices). (c) Plots of B^ (red) and By 
(blue) perpendicular to the shock front. 
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Fig. 9. (fl) Plot of Wy(0,y) for -2 < y < 2 at f = 2.6 (red) and 
t = 3.0 (blue). The perturbation has crossed the null point. The 
dashed line indicates v,,(0,3') = 0. (b) Contour of at f = 4. The 
black lines denote the separatrices. 



3.3. Oscillatory Behaviour (4.1 < f < 60) 

The evolution of for 2.9 < f < 60 can be seen in Figure |6] 
We can see that the evolution proceeds in two separate ways. 
Firstly, some of the wave now escapes the system: this propaga- 
tion can be seen above and below the location of the null point. 
Secondly, the (deformed) neutral point itself continues to col- 
lapse and forms a horizontal current sheet. This can be clearly 
seen at f = 4, and contours of are shown in Figure |9j). Here, 
current density can be seen at the four slow oblique magnetic 
shocks, along (approximately) -0.2 < x < 0.2, y = (i.e. a hor- 
izontal current sheet), at the locations of our two shock-cusps 
and due to the wave propagating away. Thus, we see that current 
density forms at several locations. The formation of a current 
sheet was not seen in [Paper I| (which reported the formation of a 
current density line at the null point). 

Returning to Figure |6] we see that this horizontal current sheet 
shortens in length (f = 4.1). This is because the jets to the left 
and right of the origin heat the plasma, which begins to expand. 
This expansion squashes the horizontal current sheet, forcing the 
separatrices apart. The (squashed) horizontal current sheet then 
returns to a 'cusp-like' null point which, due to the continuing 
expansion from the heated plasma, in turn forms a vertical cur- 
rent sheet (f = 7). The evolution now proceeds through a series 
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Fig. 10. Plot of time evolution of ;,(0,0) for < f < 60. 
Insert shows plot of time evolution of j~{Q, 0) for 25 < f < 60 
(same x-axis, different y-axis). The dashed lines indicate max- 
ima (red) and minima (blue) in the system and the green line 
shows X-(0,0) = 0.8615. 



of horizontal and vertical current sheets and displays oscillatory 
behaviour (as seen by Craig & McClymont 1991) . At f = 60, we 
see that the majority of the initial wave pulse has propagated 
away from the X-point, and that the associated velocities are 
negligible. It is also interesting to note that the final state (f = 60) 
is non-potential; the separatrices for f = are overplotted in red 
and there is clearly a small offset. This is due to the finite amount 
of current left in our system (see below). 

The oscillatory nature of the system can be clearly seen from the 
time evolution of j^{Q, 0) shown in Figure[TO] The red / blue lines 
indicate maxima / minima in the system and the green line shows 
^.(0, 0) - 0.8615, which is the limiting value of the oscillation. 
We can see that there is no current density at x = y - (i.e. 
the location of the potential null point) before f - 2.6. This is 
as expected since our initial wave pulse does not reach the null 
before this time. At f = 2.6 we see a strong spike in the current 
density: this is the result of our nonlinear wave reaching the null 
point for the first time (see Figure |5] at t - 2.6). At f = 4.4, the 
current density associated with the triangular shock-cusps (from 
the slow oblique magnetic shocks) reaches the origin; giving rise 
to a second sudden spike in current density. 
At f = 6.4, there is a third large spike in current density (opposite 
sign to previous two peaks). This indicates the formation of the 
first vertical current sheet. The peak shortly following this (at t = 
7.4) is due to the current density associated with the triangular 
shock-cusps from this vertical current sheet reaching the origin. 
The next peak in current density occurs at f = 11.8. This rep- 
resents the formation of a second horizontal current sheet (see 
Figure|6]at t = 12). However, there is no secondary peak asso- 
ciated with this horizontal current sheet, unlike the previous two 
current sheets. 

From Figure [To] we can see that at later times, the cycle of max- 
ima and minima continues, with each maxima (red lines) asso- 
ciated with the formation of vertical current sheet and each min- 
ima (blue) associated with a horizontal one. This oscillation has 
a decreasing period, and the time taken to go from a horizontal to 
vertical current sheet is shorter than vice-versa. This is because 
the system finds it easier to form vertical current sheets due to 
the asymmetric heating around the null point (see below). 
We can also see that the oscillation in current density is tending 
to a constant value of j.{Q,Q) = 0.8615. This is in agreement 
with our statement above that the final state is non-potential. This 
is because at the end of the simulation, the plasma to the left and 
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Fig. 11. (fl) Plot of P[x,Q] (red) against A,[x,Q] and plot of 
^[0,y] (blue) against Aj,[0,3']. (b) Contour of plasma pressure 
at f = 60. Contours of A-(t = 60) are overplotted in white and 
the separatrices are in black. 



right of the neutral point is hotter than above and below, due to 
the hot jets that formed and heated the (initially cold) plasma 
after f = 2. Consequently, the plasma pressure in greater to the 
left and right of the null and hence the system finds it easier 
to form vertical current sheets (due to the asymmetric heating). 
The final state (Figure|6]at t = 60) shows that the X-point is very 
slightly closed up in the vertical direction, in agreement with /- 
tending to a small positive value. 

If the final state is in force-balance, we would expect the pres- 
sure to be a function of A,, i.e. the plasma pressure should lie 
along contours of A, if P = P(Ay). Figure [TTb shows that a plot 
of A,(x,0) and A46,y) against 'P(x,0) (red) and P{x,0) (blue) 
yields a single curve, and from Figure [TTb. we can see that the 
agreement between the pressure and contours of A,(t = 60) is ex- 
cellent. This implies that although the final state is non-potential, 
it is still in force-balance. 



4. Reconnection 

From Section [331 it is clear that we have oscillatory behaviour in 
our system. However, do we have reconnection? We demonstrate 
that reconnection is occurring via the following two pieces of ev- 
idence: firstly, in Figure [12] we have plotted a selection of field- 
lines in our system at four different time slices, where red field- 
lines originate from x e [-20, -19.99], y = +20 and blue field- 
lines originate from x - +20, y e [-20,-19.99], all with equal 
separation of 0.01. The separatrices (which change in time) are 
plotted in black. In order to clearly show the results, we have 
plotted a subsection of our computational box: x,y e [-0.5, 0.5]. 
At t - 0, all the red fieldlines from the top / bottom left cor- 
ner connect to the top / bottom right corner, and all the blue 
fieldlines connect from the top left / right corner to the bot- 
tom left / right corner. The (potential) separatrices pass through 
x,y - [+20, +20] and the origin, as expected from equation 
(O, and separate the red and blue fieldlines. Thus, if we have a 
change in connectivity, we should see a red / blue fieldline cross 
the (evolving) separatrices. Since we have reflecting boundaries 
(v - 0) and the value of the vector potential. A,, on the bound- 
aries does not change with time, we can be confident that we are 
always plotting the same fieldlines in each subfigure. 
Figure[T2|(f = 4) shows our choice of magnetic fieldlines shortly 
after the formation of the first horizontal current sheet (compare 
to Figure [5] at t - 4). Here, we can see that one of our red field- 
lines has crossed the separatrices; indicating a change in connec- 
tivity. Figure [T2|(f - 8) shows our choice of magnetic fieldlines 
shortly after the formation of the first vertical current sheet, and 
now we see that some of the blue lines have crossed the separa- 
trices and the red fieldlines are once again all on the same side 
of the separatrices. At the end of our simulation (f = 60), we 
see that more blue fieldlines have crossed the separatrices. Thus, 
throughout the evolution of our system, we have several changes 
in connectivity, giving us qualitative evidence for reconnection. 
Quantitative evidence for reconnection in our system can be 
seen in Figure [13] Here, we see a plot of the time evolution of 
A-(0, 0), where changes in the vector potential at the origin in- 
dicate changes in connectivity. Figure [T3k shows a plot of the 
time evolution of A,(0, 0). Before t = 2.6, Aj.(0, 0) = as ex- 
pected (no disturbance to the null point). After this time, we see 
that A-.(0, 0) oscillates but also displays a clear trend that is tend- 
ing towards a straight fine: A, - -10"^ (8.615f + 148.434). This 
straight line is associated with the final current density remain- 
ing in our system, i.e.: 



aA- 1 , 1 
ot R,„ R„ 



(5) 



Hence, if our final state contains constant current density, say 
jc=i, then we expect A, to change linearly in time. Assuming 
j, — > jc=] and integrating equation ([5]) and comparing this to 
our straight line, we see that jc=i = 0.8615, which is in excellent 
agreement with our estimate from Figure [TOl 
Figure [T3b shows a plot of the time evolution of Ai = 
-105A,(0,0)/(8.615f+ 148.434), i.e. we have removed the 
straight line trend. We have actually plotted -A] to aid com- 
parison between Figures [T3k and [T3b . We can clearly see the 
oscillatory nature of A i (0, 0) and that it tends towards a constant 
value of -1 (which gives confidence that our straight line fit is 
appropriate). We can also see that Ai(0, 0) undergoes significant 
changes at f = 2.6,4.4 & 6.4 (overplotted in blue), where these 
times correspond to the nonlinear wave reaching the X-point, the 
first triangular shock-cusps reaching the null and the formation 
of the first vertical current sheet, respectively. 
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Fig. 12. Selection of fieldlines in our system at times f = 0, 4, 8 & 60. Red fieldlines originate from x e [-20, - 19.99], y = ±20 and 
blue fieldlines originate from x - ±20, y e [-20, 19.99], all with equal separation of 0.01. The separatrices is plotted in black. The 
same fieldines are shown in each subfigure. 
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Fig. 13. (a) Plot of the time evolution of Aj.(0,0), where we have overplotted a sti-aight line fit: A, - -10"^ (8.615f + 148.434). 
{b) Plot of the rate of change of A- at x = y = with straight line trend removed: -Ai = 10^A;(0, 0)/ (8.615f + 148.434), where 
f = 2.6,4.4 & 6.4 are shown in blue, (c) (Dashed line) Plot of the time evolution of the reconnection rate, i.e. dAJdt\ (0,0), 
(dotted line) overplot of evolution of j,(0, 0). The green line shows -/?„, 5A,/5f| (0, 0) - 0.8615. 



Figure [T3b shows us that Ai(0, 0) is continuously changing as 
it tends to -1, and thus reconnection is continuously occurring 
until this time (f ■x, 50). Thus, whenever -Ai(0, 0) oscillates 
through -1, we are increasing or decreasing flux on one side 
or the other of our final state. 

Finally, we can calculate the reconnection rate in our system, 
defined as dA~jdt\QQY and we have plotted -R„i cJAj/Sfl^yo) 
Figure [T3t . We can clearly see the reconnection rate varies 
throughout the simulation and tends towards a constant value. 
However, from equation (|5]l we see that the reconnection rate 
is the same as -jzJRm and thus is expected to tend towards 
-jc=i/Rm- Indeed we have overplotted the current density evo- 
lution (dotted fine) and the agreement is excellent. 



5. Amplitudes C=0.5 and C=2 

In this section, we examine the effect of changing the am- 
plitude, C, of our initial condition in Vj^ (equation |4|. Here, 
we consider C = 0.5 and C = 2, and the results can be 
seen in Figure [14] Comparing to Figure [T3j we can see that 
the overall behaviour is very similar that of C = 1, where 
we observe oscillatory reconnection and a cycle of horizon- 
tal and vertical current sheets. Figure [T4h shows a plot of the 
time evolution of Ac=o.5(0, 0), i.e. the flux function at the ori- 
gin for C - 0.5. We see that Ac=o.5(0, 0) oscillates but dis- 
plays a clear downwards trend that is tending to a straight line: 
A, = -10"^ (6.672f + 113.452). As explained in SectionH this 
straight line is associated with the final current density remain- 
ing in our system. Figure [T4b shows a plot of the time evolution 
of -A2 = 10^Ac=o.5(0,0)/(6.672f+ 113.452), i.e. we have de- 
trended Ac=o.5(0, 0). We can clearly see the oscillatory nature of 
Ac=o.5(0, 0). As before, the repeated changes in Ac=o.5(0, 0) im- 
plies that reconnection is continuously occurring in our system. 

We can also calculate the reconnection rate in the C = 
0.5 system, defined as dAc=o.5/dt\(Q,Q) and we have plotted 
-Rm dAc=u.5/dt\ (0,0) in Figure [T4t. We see that the reconnec- 
tion rate varies throughout the simulation and tends to a constant 
value. From equation (|5j, we see that the reconnection rate is 
also an excellent measure of the current density in the system, 
tending to a final current density of y'c=o.5 = 0.6672. Comparing 
to Figure [T3t . we see that the current evolution is comparable 
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Fig. 14. (a) Plot of the time evolution of Ac=o.5(0, 0), where we have overplotted a straight line fit: Ac=o.5 = 
-10"^(6.672f+ 113.452). (b) Plot of the rate of change of detrended Ac=o.5 x - y - 0, where -A2 - 
10^Ac=o.5(0, 0)/ (6.672f + 113.452), where t - 2.6 is shown in blue, (c) Plot of the time evolution of the reconnection rate, i.e. 
-R,„ dAc=o.5/dt\ (0, 0), where the green line indicates -R„, dAc=o.5/dt\ (0, 0) = 0.6672. (d) Plot of the time evolution of Ac=2(0, 0), 
where we have overplotted a straight line fit: Ac=2 = -10""^ (1.194f + 6.242). (e) Plot of the rate of change of detrended Ac=2 at 
X - y - 0, where -A3 - 10'*Ac=2(0, 0)/ (1.194f + 6.242), where f = 2.6 is shown in blue. (/) Plot of the time evolution of the 
reconnection rate, i.e. -R,,, dAc=2/dt\ (0,0), where the green line indicates -R„ dAc=2/dt\{Q,Q) = 1.1943. 



in nature but slightly smaller in magnitude. In addition, the final 
current left in the system is smaller than jc=i = 0.8615. 
Figure [T4h shows a plot of the time evolution of Ac=2(0, 0), i.e. 
the flux function at the origin for C = 2. As before, we see 
that Ac=2(0, 0) oscillates and displays a clear downwards trend 
fliat is tending to a sti-aight line: A, - -lO""* (1.194f + 6.242). 
Figure [T4fe shows a plot of the time evolution of -A3 - 
10'*Ac=2(0,0)/(1.194f + 6.242). Again, we can clearly see the 
oscillatory nature of the system. Figure [T4h shows the recon- 
nection rate, and therefore the current density evolution in the 
C = 2 system. Comparing to Figure [T3b. we see the current evo- 
lution is comparable in nature but slightly larger in magnitude. 
The final current left in the system, jc=2 = 1.1943, is larger than 
7c=i = 0.8615. 

We have also indicated t = 2.6 in Figures [T4b ancfl4fe. i.e. the 
time taken for the nonlinear wave to first reach the X-point in 
the C - 1 system. We see that in the C = 0.5 I C - 2 system, the 
wave reaches the X-point after / before this time, since the shock 
forms earlier and travels faster the larger the value of C. 



6. Conclusions 

This paper describes an investigation into the nature of nonlin- 
ear fast magnetoacoustic waves in the neighbourhood of a 2D 
magnetic X-point. We have solved the compressible and resis- 
tive MHD equations using a Lagrangian remap, shock capturing 
code (LARE2D). We consider a circular, sinusoidal pulse in Vj^ 
as our initial condition in velocity (equation |4]i, which naturally 
splits into two waves and we focus on the wave travelling to- 
wards the null point. Implemented damping regions remove the 
kinetic energy from the outgoing wave. Initially, we consider an 
incoming wave of amplitude C - 1. 

Between < f < 1, we find that the incoming wave propagates 
across the magnetic fieldlines and keeps its initial pulse profile 
(an annulus). The annulus contracts as the wave approaches the 
null point, and this is the same refraction behaviour reported in 
Paper I , due to the spatial variation of the (equilibrium) Alfven 
speed. We also note that the incoming wave pulse is developing 
an asymmetry: the wave peaks are propagating faster (relative to 
the footpoints) in the y-direction than in the x-direction. The 
incoming wave develops discontinuities, where in the y- / x- 
direction the wave peak / trailing footpoint is catching up with 
the leading footpoint / wave peak. The development of shocks 
around f = 1 was not reported in [Paper I| as it is a nonlinear 
effect, and arises because of our choice of a velocity initial con- 
dition: in the nonlinear regime, specifying an initial condition in 
velocity also prescribes a background velocity profile (Figure 
|4|i. This phenomenon is demonstrated for a simple system in 
Appendix |B] In addition, our initial condition in Vj^ appears to 
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excite the m - mode, but this corresponds to the m - 2 mode 
in cartesian components. 

For 1.4 < t < 2.8, we follow the asymmetry reported above in 
the development of fast oblique magnetic shock waves. By con- 
sidering the physical quantities along a line perpendicular to the 
shock front, we found an abrupt increase in density, tempera- 
ture and (consequently) pressure. It is interesting to note that the 
shock has heated the initially yS = plasma, thus creating /3 
at these locations. 

At f ss 2, we observed that the shocks above and below y - Q be- 
gan to overlap (starting at x ss +1), forming a triangular 'cusp': 
the shock-cusp. This led to the development of hot jets, which 
substantially heated the local plasma and significantly bent the 
local magnetic fieldlines. The hot jets set up slow oblique mag- 
netic shock waves emanating from the shock-cusp. In addition, 
there is evidence of slow shocks along the sides of the jet up- 
stream of the tip and we see kinks in the fieldlines at the tip of 
the jet, indicative of a fast shock. Thus, the jet heating itself is 
accomplished by a combination of slow and fast shocks. It is in- 
teresting to note that the jet has a bimodal structure consisting of 
a hot, narrow jet incased within a broader, lower temperature jet, 
which is a feature that is not predicted by steady-state reconnec- 
tion theory. 

The nonlinear wave, both the fast shocks and their overlap (the 
tails of the jets), reach the null point at f = 2.6. The shocks 
have deformed the magnetic field such that the separtrices now 
touch one another rather than intersecting at a non-zero angle 
(called 'cusp-like' by Priest & Cowley ri975l ). However, as can 
be seen at later times, the separatrices continue to evolve and 
so this osculating field structure is not sustained for any length 
of time. That the perturbation can reach the null point is a phe- 
nomenon not seen in [Paper I[ The nonlinear wave then passes 
through the null point, again entirely different behaviour to that 
seen in [Paper I| 

After t - 2.6, The evolution subsequently proceeds in two sep- 
arate ways. Firstly, some of the wave now escapes the system, 
corresponding to the wave that has passed through the null point. 
Secondly, the (deformed) null point itself continues to collapse 
and forms a horizontal cuiTent sheet. CuiTent density exists in 
this horizontal current sheet, and also at the four slow oblique 
magnetic shocks, at the location of the shock-cusps and in the 
wave propagating away (as opposed to in the linear system, 
where the current density accumulated at the null point). 
After f = 4, the evolution proceeds as follows: the jets to the 
left and right of the origin continue to heat the plasma, which in 
turn expands. This expansion squashes and shortens the horizon- 
tal current sheet, forcing the separatrices apart. The (squashed) 
horizontal current sheet then returns to a 'cusp-like' null point 
which, due to the continuing expansion from the heated plasma, 
in turn forms a vertical cuiTent sheet. The evolution then pro- 
ceeds through a series of horizontal and vertical current sheets 
and displays oscillatory behaviour (as reported by Craig & 
McClvmont lWl) . 

At the end of our simulation (f - 60), we see that the majority of 
the wave pulse has propagated away from the null point, and that 
the associated velocities are negligible. It is also interesting to 
note that our final state is non-potential: there is a finite amount 
of current density left in our system {jc=\ - 0.8615). This non- 
potential state occurs because the plasma to the left and right 
of the null point is hotter than that above and below, due to the 
hot jets that formed and heated the initially /3 - Q plasma after 
t = 2. Consequently, the plasma pressure is greater to the left and 
right and hence the final state shows that the X-point is slightly 
closed up in the vertical direction. This also explains why the 



time taken to go from a horizontal to vertical current sheet is 
shorter than the reverse: the system finds it easier to form vertical 
current sheets due to this asymmetric heating around the null. We 
also found that the pressure is approximately a function of the 
vector potential, i.e. P = PiA,), in the final state, indicating that 
although the final state is non-potential, it is still approximately 
in force balance. Of course, the system will eventually return to a 
potential state due to diff'usion, but this will occur on a far greater 
timescale than that of our simulation (fdiffusion ~ Rm = 10'^). 
We provide two pieces of evidence for reconnection in our sys- 
tem. Qualitatively, we observed changes in fieldline connectivity 
(Figure \T2\ and quantitatively we looked at the evolution of the 
vector potential at the origin. We found that A,(Q, 0) oscillates 
and displays a clear trend that was tending towards a straight 
line: A- = -10"'' (8.615f + 148.434). This straight line is associ- 
ated with the final current density left in the system. Thus, since 
we have both oscillatory behaviour in our system and evidence 
for reconnection, we conclude that the system displays oscilla- 
tory reconnection (detailed by Craig & McClymont 119911 and 
reported more recently by Murray et al. |2008l l. 
We then extended our study to look at the effect of changing the 
amplitude, C, of our initial condition (equation |4]i. Considering 
both C - 0.5 and C = 2, we found that the overall behaviour 
was very similar to the C - I study, i.e. oscillatory reconnection 
and a cycle of horizontal and vertical current sheets. Looking 
at the evolution of A,(0, 0) for the two systems, labelled Ac=o.5 
and Ac=2, we saw that both oscillated and displayed a clear 
downward tend. Again, this was associated with the final cur- 
rent left in the system, where yc=o.5 - 0.6672 was smaller than 
ic=\ - 0.8615, which were both less than jc=2 - 1.1943. Thus, 
we conclude that a larger initial amplitude results in a larger 
amount of current being left in the system at the end of the sim- 
ulation, i.e. in the non-potential final state. Note that in the linear 
limit, 7c=o.O{)i = 0, since the wave cannot reach the null point in 
a finite time. 

Thus, we can now fully answer our three original questions 
(posed in ©: 

(1) Does the fast wave now steepen to form shocks, and can 
these propagate across or escape the null? 

The nonlinear behaviour is completely different to the lin- 
ear regime. We observe the formation of both fast and slow 
oblique magnetic shocks, and the nonlinear wave can now 
cross, and thus escape, the null point. We have also seen that 
the shocks (asymmetrically) heat the plasma such that (3 i^Q. 

(2) Can the refraction effect drag enough magnetic field into the 
null to initiate X-point collapse or reconnection ? 

The nonlinear wave deforms the X-point into a 'cusp-like' 
point, which in turn collapses into a horizontal current sheet. 
Expanding plasma to the left and right squashes the horizon- 
tal current sheet and the system evolves through a series of 
horizontal and vertical current sheets. Changes in the value 
of the vector potential at the origin as well as changes in con- 
nectivity demonstrate we have reconnection in our system. 
Our final state is non-potential, but in force balance. Larger 
amplitudes in our initial condition correspond to larger val- 
ues of the final current left in the system. 

(3) Has the rate of current density accumulation changed, and 
is the null still the preferential location of wave heating? 
The current density now accumulates at many locations, such 
as along horizontal or vertical current sheets, along slow 
oblique magnetic shocks and at the location of shock-cusps. 
Current density can now also leave the system since the non- 
linear wave can escape the null point. This was not the case 
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in the linear regime, where all the current density accumu- 
lated at the null point exponentially in time. 

The aim of this paper is to contribute to the understanding of 
how nonlinear MHD waves behave in inhomogeneous, magne- 
tised environments. Future work in this area will consider the 
consequences of different initial conditions, such as utilising a 
localised increase in pressure or internal energy in a y6 + 
plasma. Such an initial condition is expected to introduce slow 
waves and create a region around the null point where the sound 
speed and Alfven speed become comparable in magnitude, and 
hence a layer where mode conversion could occur (e.g. Cally 
I200n McLaughUn & Hood 2006b ; McDougall & Hood|2007l). 
Finally, this work will also be extended to nonlinear wave be- 
haviour in the neighbourhood of a 3D null point, and compared 
to linear work by, for example, Galsgaard et al. (120031) . Pontin 
& Gal sgaard ( 120071 ), Pontin et al. ( I2007I I and McLaughUn et al. 
( I2OO8I 1. 
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Fig. A.l. Contours of for a fast wave pulse initially located at a radius r - 5, and its resultant propagation at (Alfven) times 
f = 0, 0.1, 1 & 3. At f = 0.1, we can see the initial pulse has split into two, oppositely travelling wave pulses. The black lines denote 
the separatrices and the null point is located at the origin. 



Appendix A: Linear Regime 

In this appendix, we numerically solve equations ([TJ subject to 
initial condition (|4]i and set C = 0.001. This recovers the linear 



results for the fast MA wave seen in Paper I This appendix has 
been included because (0 it will be useful to directly compare 
the linear and nonlinear systems subject to the same initial con- 
dition, and (//) Paper I described a wave pulse driven in from one 
of the boundaries, which is somewhat different from the setup 
studied in the present paper. The results described below (simu- 
lations using the LARE2D numerical code) are identical to those 
inlPaper I (simulations using a two-step Lax-WendrofF numerical 



scheme), and the results can be seen in Figure lATTI We find that 
the initial condition, as expected, splits into an outgoing wave 
and an incoming wave. We concentrate on the incoming wave. 
This wave propagates across the magnetic fieldlines and travels 
at the Alfven speed, and we identify it as a linear fast MA wave, 
in a cold plasma (J3 = 0). The linear fast wave keeps its initial 
pulse profile, i.e. an annulus, and the maximum amplitude re- 
mains at C, and the annulus contracts as the wave approaches 
the null point. Since the Alfven speed is spatially varying, a re- 
fraction effect focuses the wave into the null point. 
As the length scales decrease, there is a build up of current den- 
sity, growing exponentially in time, while the velocity remains 
finite in magnitude. Ohmic heating occurs preferentially at the 
null point. The refraction effect and the preferential heating at 
the null point are the key results for linear, yS = plasma fast 
wave propagation ( [Paper Note that since the magnetic field, 
and hence the Alfven speed, is zero at the null point, the fast 
wave cannot cross the null point and never actually reaches it. 

Appendix B: Isothermal 1D HD equations 

Consider the isothermal, one-dimensional, nonlinear hydrody- 
namic equations: 

dp 8 ^ dv dv 1 dp ^ 
— = -— (pv), — +v— = — , p-c,p . 
at ox ot ox p ox 

which is equivalent to equation ([1]) with v = vx, B = 0. Now let 
us non-dimensionalise, i.e. v = vqv*, p = pop*, P - Pop*, x - 
xqx*, t - tot*, where we let * denote a dimensionless quantity 
and vo, po, po, and to are constants with the dimensions of the 
variable they are scaling. We set vq - xo/to and Vo = Cj, i.e. non- 
dimensionalise with respect to the sound speed. For the rest of 
this appendix, we drop the star indices. Letting i/r - logp gives; 

dx 



Adding and subtracting gives: 

d d 

- (V + iA) + (V + 1) — (V + (A) = 0, 

Ot ox 

(V - ^) + (V - 1) — (V - <A) = 0. 

Ot Ox 

This can be solved in general by the Method of Characteristics: 
V = ]^F[x-{y + \)t]+^-G[x-(w-\)t], 
ifr = ^F[x-(v+l)t]-^G[x-(y-l)t], 

where F and G are arbitrary functions, specified by the initial 
conditions. 

Now let us consider an initial condition in velocity: 



( D( l+x) for -1 < .=c < 
v(x) = i D(l -x) forO < X < 1 
otherwise 



(B.l) 



where D is our initial amplitude. This has a general solution 



v(x,t)^^F[x-iv+l)t] + ^F[x-iv-l)t] 



(B.l) 



dt dx 



dv 8v dv 
dx ' dt dx 



where F(x) - v(x, 0). Thus, we can see that the general solution 
consists of two waves, moving with speeds v + 1 and v - 1, which 
can be thought of as v -H Cj and v - Cj due to our choice of non- 
dimensionalisation. We can see the evolution of initial condition 
( IB. lb in Figure IbTTI at times f = (black lines), f = 1 (blue) and 
t - 5 (red). In Figure [BTb we see that the velocity pulse (initial 
amplitude - D - 0.001) naturally splits into two waves, as ex- 
pected from the solution, each propagating in opposite directions 
with amplitude Z)/2. \{x, t) - ^F [x - {\ + 1) t] which, since v 
is small, can be thought of as v(x, f) x ^F (x - Cst) corresponds 
to a disturbance travelling in the increasing x-direction (to the 
right), whereas y(x, t) - jF [x- (v - 1) t] ^F (x + Cgt) corre- 
sponds to a disturbance travelling in the decreasing x-direction 
(to the left). This choice of D represents the linear regime. 
In Figure [BTTb we see the evolution of a velocity pulse of initial 
amplitude - D - 0.4. Again, the initial pulse splits into two op- 
positely travelling waves: one propagating to the right (located 
~ X = 5) and one to the left (located ^ x = -5). However, we 
can now clearly see the nonlinear effects: the two wave pulses 
are beginning to 'tip over' and shock, as expected for nonlinear 
waves. However, non-intuitively, the waves are both developing 
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Fig. B.l. Evolution of equation ( IB. lb at times f = (black lines), t - 1 (blue) & t - 5 (red), for three choices of initial amplitude 
(a) D - 0.001 (linear), (b) D - 0.4 & (c) Z) = -0.4. Note we have plotted -v in (c) to aid comparison. 



discontinuities, i.e. shock fronts, on the same faces. The pulse 
propagating to the right is developing a discontinuity ai x ^ 6, 
where the peak is catching up with the the leading footpoint. 
Conversely, the discontinuity developing in the pulse propagat- 
ing to the left is located at x ^ -4, i.e. where the trailing foot- 
point is catching up with the peak. 

This asymmetry in the development of discontinuities can be ex- 
plained by equation &.2\ and can be thought of as follows: when 
we set our initial condition in velocity, we are effectively pre- 
scribing a background velocity profile, and it is this profile that 
leads to asymmetry in our system (as opposed to the symmetry 
present in Figure lBTk ). A choice of D > corresponds to a pos- 
itive background velocity profile that explains the development 
of the discontinuities on the rightmost faces of both the left and 
right propagating waves. 

This explanation is confirmed by the evolution seen in Figure 
IB. lb . Here, we set D - -0.4 and have plotted -v to aid com- 
parison with {a) and (b). We see that again the initial pulse splits 
into two oppositely propagating disturbances (each with ampli- 
tude Djl) and that both these pulses are developing discontinu- 
ities on the leftmost faces. This is in agreement with the above 
interpretation, since a choice of D < corresponds to a negative 
background velocity profile. 

Interestingly, the footpoints in all three subfigures are located at 
the same locations. This is because the nonlinear effects are only 
apparent at large amplitudes and so the footpoints propagate at a 
speed of unity, i.e. Cj. 



